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Abstract 

Evolution of quadrupole deformations in sd and pf shell nuclei with mass A= 18~56 is studied 
by using deformed Skyrme Hartree-Fock (HF) model with pairing correlations. We point out that 
the quadrupole deformations of the nuclei with the isospin T=0 and T=l show strong mass number 
dependence as a clear manifestation of dynamical evolution of deformation in nuclear many-body 
systems. The competition between the deformation driving particle-vibration coupling and the 
closed shell structure is shown in a systematic study of the ratios between the proton and neutron 
deformations in nuclei with T=|T2|=1. Calculated quadrupole and hexadecapole deformations are 
compared with shell model results and available experimental data. A relation between the skin 
thickness and the intrinsic Q2 moments is also discussed. 
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I. INTRODUCTION 



A study of the mass number dependence of the deformations in nuclei is an interesting 
subject in relation to a manifestation of spontaneous symmetry breaking in quantum many- 
body systems. The effect of spontaneous symmetry breaking is a general phenomenon known 
in many fields of physics. In molecular physics, the spontaneous symmetry breaking was 
discovered by Jahn and Teller in 1937 [1]. The coupling to the quadrupole vibration is the 
main origin of the static deformation in both molecules and atomic nuclei [2]. The pairing 
correlation in atomic nuclei is known to stabilize the spherical symmetry. An unique feature 
of the evolution of deformation in atomic nuclei will appear in the competition between 
the deformation deriving particle- vibration coupling and the pairing correlations [3] . Top of 
that, the shell structure of nuclei will give rise to a variety of different shapes, i.e., prolate, 
oblate and triaxial shapes depending on the position of the Fermi energy between two closed 
shells. 

Electromagnetic observables provide useful information to study structure of nuclei of 
both ground states and excited states. Namely, the observables such as quadrupole (Q2) 
moments and electric quadrupole (E2) transitions are expected to pin down precise infor- 
mation of deformations. The evolution mechanism of deformation will change the shapes of 
nuclei from the beginning to the end of the closed shell. Up to now, the mass number depen- 
dence of the deformations has been studied mainly in rare earth nuclei. However, it is rather 
difficult in rare-earth and heavy nuclei to study the deformations systematically between two 
closed shells due to the limitation of available nuclei for experiments. On the other hand, 
it is possible to access experimentally from the beginning to the end of the closed shell in 
sd and pf shell nuclei. Thus, the sd and pf shell nuclei might be promising candidates to 
discuss the dynamical evolution of deformation since several electric quadrupole transitions 
are already observed [4] and new data of proton- rich unstable nuclei are also available [5-7] . 

The main aim of this paper is to study the dynamical evolution of nuclear deformation 
in sd and pf shell nuclei with the isospin T=0 and 1 including unstable proton-rich nu- 
clei. To this end, we perform deformed Skyrme Hartree-Fock (HF)+BCS calculations with 
a density-dependent pairing interaction in nuclei of mass number A=(16~56). Firstly, the 
mass number dependence of proton deformation is examined in comparison with empirical 
values extracted from the electric quadrupole transition strength B(E2) of the T=0 nuclei. 
The quadrupole moments of the first excited 2+ states are also studied to extract the sign 
and the magnitude of the intrinsic quadrupole deformations. Then the ratios between calcu- 
lated proton and neutron deformations of T=l nuclei are compared with shell model results 
and empirical values obtained from B(E2) values of mirror nuclei with T=l and T^ = ±1, 
assuming the isospin symmetry in the two nuclei. The evolution of higher multipole de- 
formation is also an interesting subject related with the shell filling effect of the deformed 
potential. We will study the hexadecapole deformation by using the deformed HF wave 
function and compared with available empirical data. 

This paper is organized as follows. The effect of the dynamical symmetry breaking 
on nuclear deformation is reviewed in Section llA. Deformed HF+BCS calculations with 
Skyrme forces and a density dependent pairing interactions are presented in Section IIB. 
The calculated quadrupole moments are compared with experimental data in Section IIB 
and lie The isospin symmetry of the intrinsic Q2 moments is discussed in comparison 
with shell model calculations and other mean field models in Sections 111 A and IIIB. A 
relation between the skin thickness and Q moments is mentioned in Section IIIB. Section 
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IIIC is devoted to the study of hexadecapole deformations of the sd— and p/— shell nuclei 
calculated by the HF wave functions and compared with available experimental data. A 
summary is given in Section IV. 



II. DEFORMATIONS IN MEDIUM-HEAVY NUCLEI IN MEAN FIELD THEORY 
A. Spontaneous symmetry breaking and pairing correlations 

The phenomenon of spontaneous symmetry breaking appears in many field of physics. 
As an example of the symmetry breaking, the deformation effect was observed commonly 
in quantum many-body systems, not only in molecules, but also atomic nuclei and atomic 
clusters. The spontaneous symmetry breaking was recognized also in the field theory, su- 
perconductors and condensed matter physics. Recently, it was discussed in relation to the 
chiral band of tri-axial nuclei [8] . 

Let us briefly discuss the spontaneous symmetry breaking effect to derive the nuclear 
deformation. The simplest example is the case of two degenerate single— particle states in- 
teracting with one collective mode denoted by the normal mode coordinate Q. The collective 
hamiltonian for the system can be given by the sum of the kinetic energy and the potential 



where the energy unit of this hamiltonian (1) is the eigenvalue of the normal mode. This 
collective hamiltonian is obtained based on the Born-Oppenheimer approximation, i.e., the 
decoupling between the collective and single particle motions. This hamiltonian gives the 
potential minima of = at both side of the deformation 



The positive value +k gives a prolate deformation, while the negative value —k gives an 
oblate deformation. The pairing correlations play an important role for the structure of not 
only ground states but also excited states in nuclei. The interplay between the deformation 
driving force and the pairing correlations gives an unique feature of the symmetry breaking 
of atomic nuclei compared to other quantum many-body systems [3]. We will discuss the 
effect of pairing correlations on the deformation later in this section and also in Section IIIA. 

For a harmonic oscillator potential in the limit of large quantum numbers, there is a 
symmetry between the prolate shape occurring at the beginning of the shell and the oblate 
shape at the end of the shell [9]. The HF neutron single-particle energies of ^''Si are given in 
Fig. 1 as a function of the deformation parameter /32. In sd shell nuclei, the strong preference 
for prolate deformation in the beginning of the shell can be understood from the fact that 
the energy of the orbit K'^=l/2+ from 1^5/2 orbit (in the spherical limit) decreases more 
strongly as a function of deformation than do any of the orbits in a potential of oblate shape 
near the Fermi surface as is shown in Fig. 1. Namely the strength of the linear term in 
the hamiltonian (1) is the largest for the prolate side at the beginning of the shell in the 
HF calculations. For the harmonic oscillator potential, one would expect a preference of the 
oblate shape at the end of the shell due to the symmetry of the single— particle energies. 
However, the spin— orbit force and the term in the potential reduce the driving force for 
oblate deformation although the argument of the harmonic oscillator potential prefers it in 
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the latter half region of the closed shell. Thus, in most of nuclei, the calculated and observed 
deformations are very much dominated by the prolate shape [10]. However, in light nuclei 
like C and Ne isotopes, it was pointed out that the deformation will change in a similar way 
to be expected from the deformed harmonic oscillator potential since the spin— orbit force 
and the term are relatively small in the mean fields of these isotopes [11]. 

Higher multipole deformations are also associated with the spontaneous symmetry break- 
ing of nucleus. The effect of the shell filling in deformed potentials was discussed on the 
higher multipole deformations in refs. [9] and [12]. We focus on the hexadecapole deforma- 
tion Q4 in this paper and discuss the mass number dependence of it. 

In the past, many theoretical model calculations were performed to study the masses 
and the deformations in sd and pf shell nuclei (See, for examples, review articles of various 
models in refs. [13, 14]. Ragnarsson and Sheline carried out the Nilsson Strutinsky type 
calculations with modified oscillator potential to discuss the deformation of nuclei in the 
whole nuclear periodic table [15, 16]. MoUer and collaborators performed very extensive 
calculations with a finite-range droplet model (FRDM) in the whole mass region from ^^O 
to ^^^136 [17]. Various self-consistent models are also applied to study the medium mass 
nuclei; HF model [12, 18, 19], HF+BCS model [20], HF+Bogoliubov model [21, 22]. Tri-axial 
rotor model was also used to study the deformation in [23]. 



B. Mass number dependence of deformations in T=0 nuclei 

We investigate the mass number dependence of deformations in sd and pf shell nuclei 
with the isospin T=0 and T=l. To this end, we perform deformed HF-I-BCS calculations 
with two Skyrme interactions SGH and SHI to see the interaction dependence of the results. 
Both interactions are commonly used in the mean field calculations and also random phase 
approximations for the excited states. The axial symmetry is assumed for the HF deformed 
potential. The pairing interaction is taken to be a density dependent pairing interaction in 
BCS approximation; 

V{v^,v,) = V^(^l-^^ ^{r^-r^) (3) 

where p(r) is the HF density at r = (ri-|-r2)/2 and po is chosen to be 0.16 fm~^. The pairing 
strength is taken to be Vq = —410 MeV-fm^ for both neutrons and protons [24]. A smooth 
energy cut-off is employed in the BCS calculations [25]. We study in this section quadrupole 
moments Q2 using the constrained deformed HF wave functions. The constraint is imposed 
on the axial mass quadrupole moment operator Q2o- The intrinsic proton and neutron 
Q,2n nioments are calculated by using HF wave function \K >, 

Q2p = <^l/ Pp{x,y,z){2z^-x^-y^)dr\K> (4) 

Q2n = <^l/ Pn{x,y,z){2z^-x^-y^)dr\K> (5) 

where Pp and p„ are the proton and the neutron density operators. The intrinsic proton 
hexadecapole and neutron hexadecapole Q4„ moments are also calculated by using HF 
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wave function \K 

Q4^ = <K\j p^{x,y,z){2z''-&z\x^-ry^)-r\{x^-ry^f)dv\K> (6) 

Q4n = <^\J Pn{^,y,^){^z'' - Gz\x' + y') + ^{x' + y'f)dr\K > (7) 

The electric multipole transition strength B(EA) in the laboratory system can be ex- 
pressed as 

2A + 1 

B(EX; Kh ^ Klf) = ^^(eQ;,^)^ < UK\Q\IfK >^ (8) 

with the intrinsic Qxp moment. Furthermore, the quadrupole moment of the excited 2"*" 
state Q2(2+) can be calculated by using the intrinsic quadrupole moment Q2p; 

Q2(2+) = -^Q2p- (9) 

The magnitude of intrinsic Q^p moments can be extracted experimentally from observed 
B(EA) values by using Eq. (8). The observed Q2(2"'") moment will be useful to extract not 
only the magnitude but also the sign of the intrinsic Q,2p moment from Eq. (9). 

In Fig. 2, the calculated HF energy surfaces of T= nuclei with SIII interaction are 
shown for nuclei with the mass A= 16~56. The proton, neutron and mass deformation 
parameters /?2p, P2n and P2 at the minimum of the energy surface are defined by 

Q,2n 



n 



m 



where < >p, < and < >m are the proton, neutron and mass mean square radius, 

respectively. The calculated deformation parameters are tabulated in Table I for SIII and 
SGII interactions together with the intrinsic proton and neutron quadrupole moments Q,2p 
and Q,2n- The shape isomers are also found in several nuclei at low energies, which are also 
tabulated in Table I. The sign of deformation is always the same for protons and neutrons 
in a given configuration of each nucleus, while the magnitude is largely different in several 
configurations. 

An interplay between the single particle energy gap around the Fermi surface and the 
quadrupole vibration energy plays an essential role to drive the deformation. As is expected 
from the large energy gaps at the closed shells with N=Z=8 and N=Z=20, the doubly 
closed shell nuclei ^^O and ^°Ca show spherical shapes. The prolate deformations grow in 
^°Ne and ^^Mg. In the middle of the two closed shells, the prolate and the oblate shapes are 
competing in ^^Si and ^^S so that the energy surfaces become rather flat. Just before the 
end of the closed shell N=Z=20, the oblate shape appears as the ground state configuration 
of ^^Ar. Above the mass A=40, the prolate deformations are dominant in ^^Ti, ^^Cr and 
^^Fe. Then the shape becomes again spherical in a closed shell nucleus ^^Ni. The local 
oblate minima are found at low excited energies in the three nuclei A=44,48 and 52 . Since 
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the slopes of K=l/2~ and K=3/2~ states from /7/2 orbit go downward steeper in the the 
prolate side, as can be seen in Fig. 1, the prolate deformations are favored in the ground 
states of these nuclei. On the other hand, in the nuclei A=28, 32 and 38, the whole sd shell 
orbits are involved and give the competition between oblate and prolate deformations. In 
^^Ti, it is found that the prolate and the oblate minima are almost degenerate in energy. 
However the two energy minima are very shallow and have small deformations so that the 
dynamical fluctuations will dominate low excitation energy spectra and it is unlikely to have 
any rotational bands. 

The calculated Q,2p moments are compared with experimental data in Fig. 3. The exper- 
imental Q2p moments are extracted from the transition strength B(E2) between the ground 
states and the first excited 2+ states by using Eq. (8). A clear manifestation of dynamical 
evolution of deformation can be seen in both the theoretical and experimental Q,2p moments. 
Namely, the Q,2p moment is very small in the closed shell nucleus with A=16. Then it in- 
creases when few particles are added to the closed core and shows a kind of plateau at the 
middle of the shells. A similar pattern is seen in the nuclei between A=40 and 56. Quan- 
titatively, the calculated Q,2p moments are about 10~20% smaller than the empirical one 
in nuclei with A=20— 36. However, the agreement is better in pf shell nuclei. Especially 
SGH interaction gives almost identical results to the empirical ones of nuclei with A=44— 52 
within the statistical uncertainties. 

It is shown that the calculated deformations of the ground states ^^Si and "^^Ar are oblate 
shapes. In both nuclei, the slopes of deformed HE single-particle energies in Eig. 1 favor 
the oblate deformations. As a result, in ^^Si, the oblate shape is induced by the occupation 
of K'^=l/2"'" state coming from ld5/2 state in the spherical limit, while the K'^=5/2+ state 
is going up rapidly in the prolate side. The oblate shape in ^^Ar is traced back to the 
occupation of the K''=3/2"'" state coming from ld3/2 state in the spherical limit as well. 

Deformed nuclei in the sd shell have studied with various theoretical models. As far as the 
quadrupole deformations of N=Z nuclei with A=(20~36) are concerned, the present results 
are very similar to the previous calculations by the HE models with various interactions 
[12, 13] and the modified Nilsson model [16]. Namely, ^"Nc, ^"^Mg and '^^Mg were found to be 
prolate in the ground states, while ^*Si and "^^Ar are oblate. Thus, the evolution mechanism 
of the deformation mentioned here is common in the various microscopic models. 

C. Mass number dependence of deformations in T=l nuclei 

Eig. 4 shows the binding energy surfaces of nuclei with T= 1 with T^ = —1. The proton, 
neutron and mass deformation parameters /32p, (32n and ^2 at the minimum of the energy 
surface are tabulated in Table II for SIII and SGII interactions together with the proton and 
neutron quadrupole moments Q,2p and Q,2n- To study the mirror symmetry of deformations, 
the results of nuclei with T= 1 with Tz=l are also listed in Table HI for SIII and SGII 
interactions. Since it is known that the correlations beyond the mean field approximations 
are important in nuclei near ^^O and ^''Ca [26, 27], the results of nuclei with A=18 and 
40±2 are shown to illustrate the systematic trend of the deformation from the beginning to 
the end of the shells as a manifestation of the nuclear Jahn- Teller effect. The large prolate 
deformations are found in A=22 nuclei ^^Mg and ^^Ne with (32 ~0.4. The calculated proton 
quadrupole moment Q,2p "^50 fm^ (43 fm^) for ^^Mg (^^Ne) shows good agreement with the 
empirical value Q2p= 60±11 fm^ (48.1 ±1.0 fm^ ) [4]. The calculated shapes are prolate in 
nuclei with smaller mass than the middle of the two shell closures, i.e., in ^^Mg and ^^Si. 
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The prolate and oblate minima appear at almost the same energy in ^^Si and ^^Mg, while 
the energy surface is very flat in the case of as can be seen in Fig. 4(b). In ^^Ar and '^"^S 
having heavier mass than the middle of the two closed shells, the oblate minima are found 
slightly lower than the prolate ones in energy. In the pf shell nuclei, the prolate shapes 
dominate in the nuclei with A=46, 50 and 54 as seen in Tables II, and III. 

Fig. 5 shows the ratio Q2p/Q2ri divided by the ratio of proton to neutron number Z/N 
in nuclei with T=l and T^ = —1. Because of strong proton— neutron interaction in the 
mean field, the value (Q2p/Q2n)/(Z/N) becomes close to unity when the deformation is 
well developed. The empirical neutron deformations in the nuclei with T=l T^ = — 1 are 
extracted from the proton deformations in the corresponding T=l T^^l nuclei with the 
same mass assuming a mirror symmetry of the deformations between the two nuclei. As 
far as the deformed HF calculations are concerned, the mirror symmetry of deformation 
is well conserved in the nuclei of A=(18— 54) as seen in comparisons with Tables II and 
III. The ratio (Q2p/Q2ra)/(Z/N) deviates largely from unity near the doubly closed shells 
A=16 and 40. Namely the ratio is more than 2 in A=18 system both in the HF calculations 
and the experiments. Since the driving force of deformation is only two protons in ^^Ne, a 
small proton deformation is induced. While there is no neutron particles which derive the 
neutron deformation in ^^Ne, the strong proton-neutron interaction raises smaller amount 
of deformation also for neutrons. The same trend can be seen in the nucleus with ^^Ti 
having two protons top of the doubly closed shell nucleus '^''Ca. The ratio {Q,2p/Q,2n)/ (Z/N) 
becomes opposite in the case of A=38 where two neutron holes create a driving force for 
the deformation so that the neutron deformation is larger than the proton one although the 
absolute magnitude is small. 

To conflrm theoretical conjectures of the oblate deformations, it is necessary to obtain 
empirical information of the sign of the quadrupole deformation. To this end, experimental 
data of the quadrupole moment of excited 2+ state Q2(2+) will be useful to determine not 
only the magnitude but also the sign of the intrinsic quadrupole deformation with Eq. (9). 
The present HF results of Q2(2"'') arc compared with experimental data and shell model 
calculations in Table IV. In general, the HF results show good agreement with the observed 
data in both the magnitude and the sign except the two nuclei ^^O and '^■^Ca in which the 
correlations beyond the mean field approximation are important [26]. It is seen that the 
observed data show clear evidence of the oblate deformations in ^®Si and ^^Ar. The shell 
model calculations are also given in Table IV. In sd shell nuclei, the shell model results are 
very close to the HF ones, while we notice some differences between the two models in the 
pf shell nuclei. It is remarkable that our HF results give equally good or slightly better 
results in the pf shell nuclei compared with modern version of shell model calculations. 

The two interactions give almost equivalent results in Table IV except ^°Si in which even 
the sign is different. The Q2(2"'") moments in ^°Si and ^^S were discussed in many mean field 
models [13, 18, 22] and a triaxial rotor model [23]. The HF and HF+Bogoliubov calculations 
with different interactions gave large positive (oblate) Q2(2''") in ^°Si as was the same in the 
triaxial rotor model. One exception was the SII interaction which gave negative (prolate) 
value [12]. In Table IV, the SHI gives a negative Q2(2+) which is close to the shell model 
value, while the SGII shows a small positive value. At the moment, the accuracy of the 
experimental data is not good enough to distinguish which interaction or which model is 
better or not to predict the Q2(2''') moment in ^°Si. 
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III. DISCUSSIONS 



A. Mean field models and shell models 

An extensive study of deformed HF+BCS calculations was performed in ref. [32] by 
using SIII interaction with a seniority pairing interaction. The results in ref. [32] are quite 
similar to the present ones as far as the deformations of the ground states of sd shell nuclei 
A=16~40 are concerned. On the other hand, in T=0 nuclei of p/ shell, the results of the 
two calculations are different qualitatively and also quantitatively with the same Skyrme 
interaction SIII. In ^^Ti, the present results with the surface 5 pairing interaction show a 
finite prolate deformation with /32=0.12 for SIII interaction and /92=0.19 for SGII interaction, 
while the result in ref. [32] shows no sign of deformation. For and Q2n values in ^^Cr 
and ^^Fe, the present results are (20~ 30)% larger than those of ref. [32] and close to the 
empirical values as shown in Fig. 3. For pf shell nuclei with T=l , the present results show 
finite prolate deformations (32 ~0.2 for "^^Cr and j32 ~0.12 for ^"^Ni, while the results of ref. 
[32] are spherical in both cases. The Q2 moments of the present results in ^°Fe are about 
25% larger than that of ref. [32] and show better agreement with experimental data. 

The magnitude of the quadrupole deformation depends on the strength of the pairing 
interaction. In Fig. 6, the energy surfaces of ^^Cr are shown changing the pairing strength 
multiplying a factor 0.0~1.2 by the surface type pairing interaction (3). Without pairing 
correlation, the energy minimum appears at j32 ~0.3. The energy minimum becomes shal- 
lower having smaller (52 when the pairing strength is increased. The adopted pairing strength 
(in the case of the factor 1.0) in this study gives the minimum at (32=0.2. As is expected, the 
binding energy becomes larger with stronger pairing interaction. The pairing gain energy is 
about -16MeV in ^^Cr and the calculated total binding energies with the adopted pairing (3) 
are E=-380MeV in ^^Cr which is close to the empirical one E(exp)=-382MeV. As far as the 
pairing gain energy is concerned, the pairing correlation is about 20% stronger in ref. [32] 
which gives in general smaller deformations in pf shell nuclei and predicts spherical shapes 
in some nuclei like "^^Cr and ^^Fe. 

The Skyrme HF calculations of T=0 nuclei with A=32~48 were performed in ref. [33] on 
the symmetry unrestricted basis without paring correlations. They found a similar isotope 
dependence of the quadrupole deformations to that of the present calculations with the 
pairing correlations. They pointed out also that the octupole deformation is very small in 
the ground state configurations of the T=0 nuclei. 

The proton and neutron transition matrix elements Mp, M„ were calculated by the shell 
models in sd shell nuclei by Brown and Wildenthal [28] and in pf shell nuclei by Honma et 
al.,[29]. The ratio of the Q moments (Q2p/Q2n) of the shell models is obtained by using an 
equation with the values Mp, M„ of T=l, T^ = — 1 nuclei, 

(Qp/Qn) = {e;ffM, + etffM^)/{e;ffM^ + e^^^M^), (13) 

which corresponds to the ratio of transition matrices between the mirror nuclei with T=l, 
T^ = -1 and T=l, T, = 1. The effective charges are taken to be e^ff^LSB and e^-^^=0.35 
for sd shell nuclei, and €^^^—1.5 and for pf shell nuclei. In general, the mass 

number dependence of the shell model values (Q2p/Q2n)/(Z/N) in Fig. 5 is similar to that 
of the HF results. Namely, in sd shell, the absolute values (Q2p/Q2n)/(Z/N) of the shell 
model calculations are close to those of HF calculations and show good agreement with 
empirical data except the nuclei with A=18 and 38. The shell model results deviate from 



8 



unity much larger than those of HF and also empirical ones in the nuclei near the closed 
shells N=Z=8, N=Z=20 and N=Z=28 as seen in Fig. 5. The shell model transition matrix 
elements depend largely on the adopted model space for nuclei near the closed shell. In the 
case of A=18, a complete basis of (lpi/2, lc?i/2, 2si/2) model space (so called ZBM basis) [34] 
gives a better result for the value (Q2p/Q2n)/(Z/N) to be 1.67, instead of more than 3 in 
ref. [28], in comparison with the empirical data. 

In pf shell nuclei, the shell model results are close to those of HF calculations and 
consistent with the empirical values in the nuclei with A=46 and 50. On the other hand, 
the shell model values are much larger in A=42 and smaller in A=54 in comparison with 
the HF results. These results suggest that a strong proton— neutron interaction in the HF 
calculations induces smaller asymmetries in the deformations between T=l, T^ = ±1 nuclei 
near the closed shells than expected from the shell model results. In the middle of the 
closed shells, both the HF and shell model calculations give reasonable collective transition 
strength compared with the experimental values. 

Finite-range droplet model (FRDM) has been known to provide useful information on 
the masses and the deformations of nuclei in a very wide region of the mass table [17]. 
Compared with the HF model, the FRDM gives quite different results for the quadrupole 
deformations in many sd and pf shell nuclei. For examples, the FRDM gives large defor- 
mations /92=(0.32~0.37) for ^°'^^Ne and ^^'^^Mg which arc consistent with the HF results in 
Tables I, II, III. A nucleus ^^Mg is predicted as an oblate shape in FRDM, while the oblate 
and the prolate minima are almost degenerate in the HF results in Table HI. The S, Ar, 
Cr and Fe isotopes show substantial deformations in the HF calculations, while these iso- 
topes are predicted to be spherical in the FRDM. It seems that the FRDM might not give 
proper deformations in these medium mass nuclei in comparison with empirical information 
discussed in this paper. 

B. New Experimental Data and Comparisons with Theories 

Recently, the B(E2) values of several proton-rich unstable pf shell nuclei with T=l and 

T2=— 1 are observed by Coulomb excitations [5-7]. The mirror symmetry in T=l nuclei 
with A=46,50 and 54 is experimentally confirmed by these new data. Namely, in A=46, 
the experimental values are B(E2) = (950±50) e^fm"^ in ^^Ti and (929 ±199) eVm^ in ^^Cr, 
while the calculated values are B(E2)=582.1eV^^ and 913.7 e'^fm'^ for ^^Ti and ^^Cr, 
respectively, with SGII interaction. The different calculated Q moments of the two nuclei 
are due to the proton skin effect in "^^Cr which enhances the value as is expected from 
Eq. 10. Experimental data of T=l and Tz=—1 nuclei with A=50 and 54 were reported to 
be B(E2) = (1359 ±261) e^fm"^ for ^^Fe and B(E2) = (590±168) e^fm"^ for ^'^Ni in ref. [7]. 
Corresponding experimental values for T=l and T^^l nuclei are B(E2)=(1080 ±60) fm^ 
for 5°Cr and B(E2) = (620 ±50) e^fm^ for ^^Fe [4], respectively. Calculated HF values of 
A=50 nuclei with SGII interaction, B(E2)=1160e2/m^ for ^°Fe and 1097e^ fm'^ for ^oCr, 
show a clear mirror symmetry and give good accounts of the experimental data. The mirror 
symmetry is also seen in the HF calculations of A=54 nuclei with B(E2)=373e^/m^ for ^^Ni 
and 390e^/m^ for ^^Fe, although the absolute values are smaller than the observed ones. In 
refs. [5, 6], the B(E2) values for the ground states to the first 2+ states were observed by 
the intermediate Coulomb excitations which give almost identical results given in refs. [4] 
and [7]; B(E2) = (626±169) e^fm"^ for ^^Ni and (640 ±13) e^fm"^ for ^^Pe. 

Shell model calculations of B(E2) strength from the ground states to the first excited 
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states of A=50 and 54 nuclei were performed in a restricted pf shell model space with 
GXPFl effective interaction [29]. Taking the effective charges ep=1.5 and e„=0.5, the shell 
model calculations give a good mirror symmetry in the case of A=50 nuclei with B(E2)= 910 
e^/m^ for ^'^Cr and B(E2)=954e^/m^ for ^°Fe although the absolute values are somewhat 
smaller than the experimental values and those of deformed HF calculations with SGII 
interaction. The shell model gives a reasonable B(E2) value for ^^Fe to be B(E2)=651e^/m^ 
, while the calculated value is smaller than the experimental one in the case of ^^Ni as 
B(E2)=324eVrn^ 

If the mean square radii are close < >p—< >n in a nucleus, the ratio Q2p/Q2n 
divided by Z/N should be the same as the ratio Sp/Sn (or equivalently One can 

notice that T=0 nuclei in Table I show the proportionality between Q moments and the 
deformation parameters (32 as is expected from Eqs. 10 and 11. On the other hand, in 
Tables II and III, the Q moments are not always proportional to the Pp and /3„ values. This 
is due to a difference of the proton skin effect in the T^ = — 1 nuclei to the neutron skin 
effect in the T^ = 1 nuclei. In Fig. 8, the calculated proton skin thickness of the ground 
state, 

SVpn =< Tp > - < r„ >, (14) 

is shown for the nuclei in sd and pf shell with T=l and T^ = ±1. The results of SGII are 
shown in Fig. 8. We calculated also the skin thickness by using the SIII interaction and 
found almost identical results to those shown in Fig. 8. In A=18, the proton-rich nucleus 
^^Ne has a large proton skin and the neutron-rich nucleus ^^O has an equally large neutron 
skin. The large skin thickness is created by two particle occupation of the 1(1^/2 orbit which 
has much larger radius than Ip orbits. The proton skin thickness of nuclei with T^ = — 1 
decreases gradually for heavier pf shell nuclei, but still finite to be Svpn ~ 0.1 even in ^^Ni. 
The irregularity at A=30 is due to the occupation of 2si/2 orbit. The neutron skin thickness 
of nuclei with T^ =1 decreases more rapidly for heavier sd shell nuclei and almost disappears 
in the pf shell nuclei. This difference between the nuclei with T^ = ±1 is due to the effect 
of the Coulomb interaction on the mean field proton potential. Since the Coulomb potential 
makes shallower proton HF potential in pf shell nuclei, the extra two protons give always 
the enhancement of the proton rms radii in T^ = — 1 nuclei. On the other hand, the neutron 
potential is deeper than the proton one in T^ =1 nuclei so that the extra two neutrons do 
not create appreciable neutron skin in pf shell nuclei. This difference of the skin thickness 
gives rise to a substantial effect on the Q moment. One of the clear examples of the skin 
effect is seen in ^°Fe where the proton deformation is smaller the neutron one, but the Q,2p 
moment is larger than moment. 

C. Hexadecapole deformation in sd and pf shell nuclei 

Higher multipole deformations such as the hexadecapole deformation were discussed in 
ref. [12] by the Skyrme HF calculations without pairing correlations. It is an interesting 
subject whether the pairing correlations are important or not for the evolution of higher 
multipole deformations than the quadrupole. The proton and neutron hexadecapole defor- 
mations of T—0 and T—1 nuclei are listed in Tables I, II and III. In sd shell nuclei, a large 
Q4 in ^°Ne and a small Q4 in ^^Mg are expected as the shell fiUing effect in prolate deformed 
potential [9, 12]. The present HF results are consistent with the shell filling effect. The two 
interactions SIII and SGII give a large difference in the Q4 moment of ^*Si. In comparison 
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with experimental data in Fig. 7, SGII gives a better prediction in ^^Si. In ^^S and ^^Ar, the 
calculated Q4 moments are relatively small compared with other sd shell nuclei, while the 
shell model gives several times larger Q4 moments in these two nuclei [37]. We need more 
data, for example, by electron scatterings, to compare with these calculated results in ^^S 
and ^^Ar, although the proton scattering data of ^^S show an indication of large value 
[38]. 

For pf shell nuclei, the Q4 moments increase drastically from ^^Ti and peaked at ^®Cr. 
In the heavier pf shell nuclei, we can see a large drop in magnitude of Q4 moment from 
^^Cr to ^^Ni. In general, SGII gives larger Q4 moments in pf shell than SIII. In sd shell 
nuclei, the HF results underestimate somewhat the empirical values, while the calculations 
give reasonable values of Q4 moment in pf shell nuclei, ^^Ti and ^''Cr. 

IV. SUMMARY 

We pointed out the strong mass number dependence of quadrupole deformations in the 
nuclei with mass A=(16~56) as a clear manifestation of the evolution of nuclear deformation 
in atomic nuclei by using the deformed HF+BCS calculations. The effect in the medium- 
heavy nuclei is unique to compare with that in rare-earth nuclei since the prolate and oblate 
deformations appear clearly in the ground states depending on the Fermi energy in the 
deformed single-particle energy levels. It is shown that the deformed HF-I-BCS model is 
successful to describe observed Q2 and Q4 moments of the nuclei with the isospin T=0 and 
T=l in sd and pf shell in the same level as modern shell model calculations. The isospin 
symmetry of quadrupole deformations of mirror nuclei with T=l, T^ = ±1 is studied in 
comparison with the empirical data and the shell model calculations. It was shown that 
the value (Q2p/Q2n)/(Z/N) is close to unity in the well-deformed nuclei in the middle of the 
shell, while the values deviate largely from the unity in nuclei near the closed shell both in 
the HF calculations and also in the empirical values. The HF results show oblate shapes in 
some sd shell nuclei and the available observed Q2(2+) moments of the first excited states 
are consistent with the theoretical predictions. The effect of the pairing interaction on the 
deformation is carefully examined by changing the pairing strength. We point out also the 
effect of large proton skin thickness on the Q2 moments in the N~Z nuclei. 
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Figures 




FIG. 1: Neutron HF single particle energies as a function of the deformation parameter /?2 
The deformed HF calculations are performed by using SIII interaction. 
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FIG. 2: Deformed HF + BCS calculations of T=0 nuclei with SIII interaction; (a) for ^^O and 

20Ne, (b) for ^4 Mg and ^SSi , (c) for and ^^Ar, (d) for ^SQa and ^^Ti, (e) for ^^Cr and ^^Fe 
and (f) for ^^Ni. The density dependent pairing interaction (3) is adopted in the calculations. 
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FIG. 3: Absolute values of proton quadrupole moments in sd and pf shell nuclei with N=Z. 
Deformed HF + BCS calculations are performed with SIII and SGII interactions. The density 
dependent pairing interaction (3) is adopted in the calculations. The experimental data are taken 
from refs. [4, 7]. 
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FIG. 4: Deformed HF + BCS calculations of T=l with SIII interaction; (a) for ^^Ne and ^^Mg, 
(b) for 26si and 3°S, (c) for ^^Ar and ^^Ca, (d) for ^^Ti and ^^Cr, (e) for ^Ope and (f) for ^^Wi. The 
density dependent pairing interaction (3) is adopted in the calculations. 
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FIG. 5: Ratio of proton to neutron quadrupole moments (Q2p/Q2n) divided by that of proton 
to neutron numbers (Z/N) in sd and pf shell nuclei with T=l and = —1. Deformed HF + 
BCS calculations are performed with SIII and SGII interactions. The density dependent pairing 
interaction (3) is adopted in the HF calculations. The shell model values are calculated by using 
Eq. (13). The shell model transition matrices of sd shell are taken from ref. [28], while those of pf 
shell are taken from ref. [29]. The empirical values are obtained by assuming a mirror symmetry 
between proton and neutron quadrupole moments in T=l and T^ = ±1 nuclei with the same mass 
number A. Experimental data are taken from refs. [4, 7] . See the text for details. 
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FIG. 6: Deformed HF + BCS calculations of °Cr with different pairing strength multiplying a 
factor 0.0~1.2 to the surface type pairing interaction (3). The SIII interaction is used. 
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FIG. 7: Proton hexadecapole moment Q^p defined by Eq. (6). The SIII and SGII interactions are 
used in the deformed HF calculations with a surface pairing interaction (3). Experimental data 
are taken from refs. [35, 36]. See the text for details. 
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FIG. 8: The difference between proton and neutron rms radii defined by Eq. (14). The SGll 
interaction is used in the deformed HF calculations with a surface pairing interaction (3) . See the 
text for details. 
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TABLE I: Quadrupole deformation parameters P2, quadrupole moments Q2 and hexadecapole mo- 
ments Q4 in sd and pf shell nuclei with the isospin T=0. The deformed HF+BCS calculations are 
performed with SIII and SGII interactions together with the density dependent pairing interaction 
(3). The axial symmetry is assumed in the HF calculations. The proton, neutron and mass defor- 
mation parameters (P2p, P2p and P2) are obtained at the energy minima of protons, neutron and 
mass potentials, respectively. The proton and neutron quadrupole and hexadecapole moments are 
calculated by using deformed HF wave functions at each energy minima. The local minima of the 
energy surfaces at low excitation energies arc also listed. 



CiTTT 
oiii 


J\ 


Uvuergy 


P2p 


P2n 


P2 


On 


^2n 










I ivie V ) 








[jm ) 


(jm ) 




[jm ) 


i6r) 


n+ 
u 


n n 


n n 
u.u 


n n 
u.u 


n n 

U.U 


n ni 

-U.Ul 


n m 

-U.Ul 


n 1 
-u.io 


n 1 

-U. 10 


iv e 


u 


n n 
u.u 


U.oo 


n Q/i 


n Qfi 
u.oo 


oy.yo 


87 
00.0 ( 


zo / .00 


9/1 Q R'i 
Z4o.0t3 




n+ 
u 


i.yo 


-U. iO 


n 1 ^ 

-U. iO 


n 1 ^ 

-U. iO 


1 9 Q/l 

- i z . y4 


1 9 
-iZ.OO 


91 '?9 
Zi.OZ 


9n n^ 

ZU.Uo 


IVlg 


n+ 
u 


n n 
u.u 


n /in 

U.4U 


n /in 

U.4U 


n /in 

U.4U 


00. ( i 


04. 4o 


1 9fi R"? 
iZD.DO 


1 9Q 1^7 
IZo.O < 


28 o,- 


n+ 
u 


n n 
u.u 


-U.Z4 


n 9/1 

-U.Z4 


n 9/1 

-U.Z4 


-oy.y / 


-00.00 


1 97 R1 
iZ ( .Di 


1 9n 

iZU. ( 




n+ 

U 


u.ou 


U.U4: 


U.U4 


n OA 

U.U4 


1 n 

U. lU 


^ Q9 


7 98 
/ .zo 


U. / 4 


32s 


0+ 


0.0 


0.20 


0.20 


0.20 


41.62 


40.10 


-49.84 


-44.64 




0+ 


1.54 


-0.08 


-0.08 


-0.08 


-16.40 


-15.59 


-8.82 


-7.84 




0+ 


0.0 


-0.18 


-0.18 


-0.18 


-44.46 


-43.01 


-47.34 


-41.74 




0+ 


2.30 


0.09 


0.09 


0.09 


22.59 


22.01 


-15.62 


-13.92 




0+ 


0.0 


0.01 


0.01 


0.01 


2.72 


2.64 


1.17 


0.97 




0+ 


0.0 


0.15 


0.15 


0.15 


51.73 


49.96 


508.84 


467.90 




0+ 


0.20 


-0.11 


-0.11 


-0.11 


-35.78 


-34.80 


155.90 


144.70 




0+ 


0.0 


0.26 


0.26 


0.26 


101.80 


98.72 


785.70 


737.63 




0+ 


3.55 


-0.11 


-0.11 


-0.11 


-41.17 


-39.86 


40.45 


35.71 




0+ 


0.0 


0.20 


0.20 


0.20 


90.72 


88.19 


-27.00 


-16.71 




0+ 


3.04 


-0.08 


-0.08 


-0.08 


-34.16 


-33.22 


-32.77 


-29.07 


56 


0+ 


0.0 


0.01 


0.01 


0.01 


4.78 


4.63 


-12.13 


-11.04 



22 





7>^7r 
A 


Energy 


Q 

P2p 


Q 

P2n 


P2 


r\ 

H2p 


r\ 


W4p 


V4n 






(^IVieV j 








( /^2\ 

[Jm ) 


( f^2\ 


[Jm ) 


[Jm ) 


U 




u.u 


u.u 


u.u 


0.0 


-0.01 


n m 

-0.01 


n 1 '7 

-0.17 


-0.17 


20 AT 




U.U 


U.41 


U.41 


u.41 


/I o o o 

43. J8 


42.12 


on'? nn 

297.09 


279.87 




n+ 
U 


z.by 


-U.iO 




-U.lb 


1 /I CI 

-14. oi 


-14.00 


on 1 
z9.Di 


0*7 /?n 

z / .bU 


24 71 ^- 

Mg 




u.u 


n /I o 

U.4z 


U.42 


U.4z 


00.05 


55.42 


129.01 


1 O A A 

125.44 


28 o- 




U.U 


r\ on 

-U./9 


-U.29 


-U.z9 


-49. Iz 


-47.58 


1 nn no 

19U.U8 


1 'TO C O 

178.53 




U 


1 n ~ 

i.;J / 


n in 
U.iU 


CI 1 n 
U.IU 


U.iU 




"1 1 f: "7 
i 1.0 / 


zi. iU 


2U.U1 


32 c 


n+ 


n n 

u.u 


U.i4 


U.io 


U.13 


27.04 


25. o3 


OT 1 

-27.16 


-24.14 


36^ J, 


n+ 


u.u 


-U.17 


-U.17 


-0.17 


-42.21 


/in /^/^ 

-40.66 


-64.86 


-57.59 






1.00 


u.uy 


U.U9 


U.U9 


22.45 


o 1 or 

21.85 


-15.14 


1 O O/^ 

-13.26 




n+ 


u.u 


u.u 


u.u 


0.0 


r\ no 

0.03 


n n /I 

0.04 


n n CT 

U.U5 


n no 

0.03 


44a^ ' 


n+ 

0^ 


u.u 


U.19 


0.19 


0.19 


65.48 


62.80 


823. Ul 


750.34 






U.DO 


-U.IU 


-U.IU 


-0.10 


-32.54 


o 1 r n 

-31.59 


14U.74 


1 on 1 

130.16 


48/^^ 


n+ 
U 


U.U 


U.oi 


U.ol 


U.ol 


123.69 


119.10 


105 / .68 


no o 1 
985.81 




0+ 


4.46 


-U.12 


-0.12 


-0.12 


-44.72 


-43.38 


52.80 


50.20 


52^g 


0+ 


U.U 


U.23 


0.23 


0.23 


99.0 


95.85 


49.32 


55.63 




0+ 


3.72 


-U.12 


-0.12 


-0.12 


-51.23 


-49.57 


106.35 


98.11 


56 


0+ 


U.U 


0.0 


U.U 


0.0 


-0.01 


-0.0 


-60.41 


-54.23 



23 



TABLE II: Quadrupole deformation parameters P2, quadrupole moments Q2 and hexadecapole 
moments Q4 in sd and pf shell nuclei with the isospin T=l and = —1. The deformed HF+BCS 
calculations arc performed with SIII and SGII interactions together with the density-dependent 
pairing interaction. See the caption to Table I for details. 
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TABLE III: Quadrupole deformation parameters P2, quadrupole moments Q2 and hexadecapole 
moments Q4 in sd and pf shell nuclei with the isospin T=l and = 1. The deformed HF+BCS 
calculations arc performed with SIII and SGII interactions together with the density-dependent 
pairing interaction. See the caption to Table I for details. 



SIII 




Energy 


P2p 




h 


Q2p 


Q,2n 


Clip 


Q4n 






(MeV) 
















ISq 


0+ 


0.0 


-0.01 


-0.01 


-0.01 


-0.51 


-1.00 


-0.13 


-0.10 




0+ 


0.0 


0.39 


0.38 


0.39 


41.38 


51.15 


206.75 


159.93 




0+ 


3.08 


-0.16 


-0.15 


-0.16 


-15.99 


-19.55 


30.64 


30.17 




0+ 


0.0 


0.29 


0.22 


0.25 


39.64 


35.41 


59.63 


64.73 




0+ 


0.31 


-0.21 


-0.21 


-0.21 


-28.33 


-33.78 


66.75 


84.20 




0+ 


0.0 


0.08 


0.09 


0.08 


13.99 


16.55 


9.92 


3.22 




0+ 


0.0 


-0.16 


-0.16 


-0.16 


-33.34 


-38.41 


-13.58 


-29.72 




0+ 


0.29 


0.16 


0.12 


0.14 


32.59 


28.67 


-27.90 


-24.61 


38^^ 


0+ 


0.0 


-0.01 


-0.01 


-0.01 


-3.49 


-1.78 


-0.39 


-0.44 


AO r^i 


0+ 


0.0 


-0.01 


-0.01 


-0.01 


-2.19 


-3.85 


-0.88 


-2.47 




0+ 


0.0 


0.20 


0.22 


0.21 


70.40 


84.65 


689.82 


688.03 




0+ 


1.29 


-0.11 


-0.11 


-0.11 


-36.09 


-39.59 


121.36 


86.17 




0+ 


0.0 


0.24 


0.23 


0.23 


93.70 


97.15 


461.74 


302.57 




0+ 


2.81 


-0.10 


-0.10 


-0.11 


-42.66 


-43.97 


43.84 


46.20 


54^g 


0+ 


0.0 


0.12 


0.09 


0.11 


52.72 


43.51 


17.18 


59.43 




0+ 


0.0 


0.11 


0.14 


0.13 


56.09 


74.64 


238.35 


343.72 




0+ 


0.17 


-0.09 


-0.10 


-0.09 


-42.38 


-50.42 


115.25 


94.07 


SGII 




Energy 










Q2n 




0,4:71 






(MeV) 
















18q 


0+ 


0.0 


0.01 


0.01 


0.01 


0.37 


1.27 


-0.19 


-0.27 




0+ 


0.0 


0.41 


0.40 


0.41 


43.19 


53.21 


225.14 


171.92 




0+ 


3.61 


-0.18 


-0.19 


-0.18 


-18.39 


-23.47 


43.48 


47.64 




0+ 


0.0 


0.32 


0.24 


0.28 


42.59 


37.48 


71.51 


73.05 




0+ 


0.16 


-0.24 


-0.24 


-0.24 


-31.71 


-38.40 


89.68 


109.92 




0+ 


0.0 


-0.02 


-0.02 


-0.02 


-3.44 


-3.57 


-5.25 


-6.15 


345 


0+ 


0.0 


-0.13 


-0.14 


-0.14 


-26.93 


-32.46 


-27.97 


-36.41 




0+ 


0.25 


0.14 


0.11 


0.12 


27.51 


25.70 


-26.42 


-19.14 


38^^ 


0+ 


0.0 


0.01 


0.01 


0.01 


2.69 


2.11 


0.38 


0.15 




0+ 


0.0 


0.01 


0.01 


0.01 


2.0 


3.83 


0.64 


1.78 




0+ 


0.0 


0.22 


0.25 


0.24 


76.50 


92.14 


807.53 


796.05 




0+ 


2.07 


-0.11 


-0.11 


-0.11 


-35.86 


-39.73 


139.57 


105.34 


50 Cr 


0+ 


0.0 


0.27 


0.26 


0.26 


105.20 


108.17 


609.74 


419.16 




0+ 


3.65 


-0.14 


-0.14 


-0.14 


-53.17 


-56.62 


94.77 


123.07 


54^g 


0+ 


0.0 


0.15 


0.11 


0.13 


62.84 


50.50 


23.90 


82.27 




0+ 


0.0 


-0.11 


-0.12 


-0.11 


-53.14 


-59.25 


173.44 


118.77 




0+ 


0.02 


0.12 


0.14 


0.13 


56.64 


73.71 


195.14 


238.98 
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TABLE IV: Quadrupole moments Q2(2^) of the first excited 2+ states in sd and pf shell nuclei. 
The unit is e • fm?. Shell model values of sd shell nuclei are taken from ref. [30], while those of pf 
shell nuclei are taken from rcf. [29]. Expcriuwnilal data arc lakcii from the compilation of ref. [31]. 



nucleus 


SIII 


SGII 


shell model 


Expt. 


180 


0.15 


-0.10 




-2.0 


-3.9 ± 0.9 


20Ne 


-11.4 


-12.4 




-12.1 


-23 ± 3 


22Ne 


-12.3 


-12.3 




-13.6 


-19 ± 4 


o 


-15.9 


-16.3 




-15.0 


-18 ±2 


28Si 


11.4 


14.0 




14.3 


16± 3 


30 Si 


-4 n 


n 97 




-6.6 


-1- 6 


32s 


-11.9 


-7.7 




-13.6 


-15.4 ± 2.0 


34s 


9.5 


7.7 




6.7 


4 ± 3 


36 Ar 


13.4 


12.1 




14.3 


11 ± 6 


42Ca 


0.63 


0.57 




1.90 


-19± 8 


46^1 


-20.1 


-21.9 




-11.1 


-21±6 


50 cr 


-26.7 


-30.0 




-26.4 


-36 ±7 


s^Fe 


-15.1 


-17.9 




-22.6 


-5 ±14 


58Ni 


-16.0 


-16.2 




-2.4 


-10 ± 6 
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